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1. Introduction 

Jet production observables are studied extensively at hadron colliders. Since the distribu- 
tion of final state jets relates directly to the parton-level dynamics, jet observables can be 
used for precision studies of QCD [1], especially in view of determinations of the strong 
coupling constant and the parton distribution functions in the proton. Experimental mea- 
surements of these observables at the Tevatron [2,3] attained an accuracy of a few per cent 
(or even better in certain kinematical ranges), and first results from the LHC [4,5] already 
show the potential for precision jet physics. Consequently, meaningful precision studies 
must rely on theoretical predictions accurate to the same level. In perturbative QCD, this 
precision usually requires corrections at next-to-next-to-leading order (NNLO). 

NNLO calculations of observables with n jets in the final state require several ingre- 
dients: the two-loop corrected n-parton matrix elements, the one-loop corrected (n + 1)- 
parton matrix elements, and the tree-level (n + 2)-parton matrix elements. For most 
massless jet observables of phenomenological interest, these matrix elements are available 
for some time already. 

The (n + l)-parton and (n + 2)-parton matrix elements contribute to n jet observables 
at NNLO if the extra partons are unresolved or are clustered to form an n-jet final state. 
Consequently, these extra partons are unconstrained in the soft and collinear regions, and 
yield infrared divergences. In these cases, the infrared singular parts of the matrix elements 
need to be extracted and integrated over the phase space appropriate to the unresolved 
configuration to make the infrared pole structure explicit. The single soft and collinear 
limits of one-loop matrix elements [6-10] and the double unresolved limits of tree-level 
matrix elements [11-14] are process-independent, and result in a factorization into an 
unresolved factor times a matrix element of lower multiplicity. 



To determine the contribution to NNLO jet observables from these configurations, 
one has to find subtraction terms which coincide with the full matrix element and are 
still sufficiently simple to be integrated analytically in order to cancel their infrared pole 
structure with the two-loop virtual contribution. Often starting from systematic methods 
for subtraction at NLO [15-18], several NNLO subtraction methods have been proposed 
in the literature [19-26], and are worked out to a varying level of sophistication. 

For observables with partons only in the final state, an NNLO subtraction formalism, 
antenna subtraction, has been derived in [27]. The antenna subtraction formalism con- 
structs the subtraction terms from antenna functions. Each antenna function encapsulates 
all singular limits due to the emission of one or two unresolved partons between two colour- 
connected hard radiator partons. This construction exploits the universal factorization of 
matrix elements and phase space in all unresolved limits. The antenna functions are de- 
rived systematically from physical matrix elements [28]. This formalism has been applied 
in the derivation of NNLO corrections to three-jet production in electron-positron annihi- 
lation [29,30] and related event shapes [31,32], which were used subsequently in precision 
determinations of the strong coupling constant [33-38]. The formalism can be extended 
to include parton showers at higher orders [39], thereby offering a process-independent 
matching of fixed-order calculations and logarithmic resummations [34,35,40,41], which is 
done on a case-by-case basis for individual observables [42] up to now. The formalism can 
be extended to include massive fermions [43]. 

For processes with initial-state partons, antenna subtraction has been fully worked out 
only to NLO so far [44]. In this case, one encounters two new types of antenna functions, 
initial-final antenna functions with one radiator parton in the initial state, and initial-initial 
antenna functions with both radiator partons in the initial state. The framework for the 
construction of NNLO antenna subtraction terms involving one or two partons in the initial 
state has been set up in [45] in the context of a proof-of-principle implementation of the 
contribution of the gg — ?• ig tree-level subprocess to di-jet production at hadron colliders. 
The initial-final and initial-initial antenna functions appearing in the NNLO subtraction 
terms are obtained from crossing the final-final antennae. Their integration has to be 
performed over the appropriate phase space. In the case of the initial-final antennae, this 
has been accomplished in [46] . For the initial- initial tree- level double real radiation antenna 
functions, partial results have been obtained in [47]. It is the aim of the present paper to 
derive the setup for NNLO antenna subtraction for single unresolved singularities at one- 
loop and to compute the integrated one-loop initial-initial antenna functions required in 
this context. 

Other approaches to perform NNLO calculations of exclusive observables with initial 
state partons are the use of sector decomposition and a subtraction method based on the 
transverse momentum structure of the final state. The sector decomposition algorithm [48] 
analytically decomposes both phase space and loop integrals into their Laurent expansion 
in dimensional regularization, and performs a subsequent numerical computation of the 
coefficients of this expansion. Using this formalism, NNLO results were obtained for Higgs 
production [49] and vector boson production [50] at hadron colliders. Both reactions were 
equally computed independently [51] using an NNLO subtraction formalism exploiting the 



specific transverse momentum structure of these observables [24], wliich was also applied 
most recently to compute NNLO corrections to associated WH production [52]. 

This paper is structured as follows: in Section 2, we construct the subtraction terms 
required at NNLO for initial-initial configurations with one unresolved parton at one loop. 
They require one-loop 2 — )■ 2 antenna functions with two partons in the initial state and 
one parton and one ofF-shell neutral current in the final state. The analytic integration of 
the initial-initial one-loop antenna functions is described in Section 3. Finally, we conclude 
with an outlook in Section 4. 

2. Initial-intial antenna subtraction at NNLO 

Antenna subtraction of initial-initial configurations at NLO is derived in detail in [44]. 
Subtraction terms with two hard partons in the initial state are built along the same lines 
as in the final-final and initial-final case. The NLO antenna subtraction term for an m-jet 
production process, to be convoluted with the appropriate parton distribution functions 
for the initial state partons, for a configuration with the two hard emitters in the initial 
state (partons i and k with momenta pi and P2), reads: 

cia^^ii) = 7y^y~^ d$m+i(/ci,... ,kj^i,kj,kj+i,...,km+i;pi,P2) -^ 

rn+l 

2 

y^ X^f.j{pi,P2, kj) Mm(ki, ■■■ , kj^i,kj^i, ... , km+i;xipi, X2P2) 
3 

Xj^ (fci, . . . , Kj_l, fcj + l, . . . , fcj„_|_l j . (,2.ij 



All the momenta in the arguments of the reduced matrix elements and the jet functions 
are redefined, which is a consequence of requiring the correct collinear factorization prop- 
erties in both initial-state collinear limits. Its should be noted that the jet function Jm 
(constructing m jets from m partons) requires all redefined momenta to be resolved. The 
two hard radiators are simply rescaled by factors x\ and X2 respectively. The spectator 
momenta are boosted by a Lorentz transformation onto the new set of momenta {/c;, / 7^ j}. 
The mapping must be based on a factorization of the (m -|- l)-particle phase space, must 
satisfy overall momentum conservation and keep the mapped momenta on the mass shell. 
In this case, this turns out to severely restrict the possible mappings. 

The tree- level antenna function X^^ ■ depends only on the incoming momenta pi,P2 
and on the outgoing momentum kj . It accounts for all singular configurations where parton 
j is unresolved and colour-connected to partons i (incoming with momentum pi) and k 
(incoming with momentum P2). The jet function Jm and the reduced matrix element in 
(2.1) depend only on the redefined momenta. With a suitable factorization of the phase 
space [44], one can perform the integration of the antenna function analytically. 

The factorization of the phase space is obtained by requiring that the two mapped 
initial state momenta should be of the form 

Pi = xipi P2 = X2P2 , (2.2) 



so that 

is in the beam axis. Since the vector component of q = pi+p2 — kj is in general not along 
the pi — p2 axis we need to boost all the momenta so that q = Aq and ki = Afc/ in order to 
restore momentum conservation. By requiring this boost to be only transverse, the phase 
space mapping is determined uniquely, resulting in the factorization 

d^rn+likl,- . . , km+i;Pl,P2) = d^rnih, ■■■ , kj-l, kj+i, ..., km+r, XiPi, X2P2) 

x5{xi — xi) 6{x2 — X2) [dkj] dxi d2;2 . (2-3) 

where [dk] = d'^k6^+\k'^)/{2TTY-'^ and 



_ Si2 - Sj2 Si2 - Sij - Sj2 
Sl2 Si2 - Sij 

X2 = ('J^^Ill 211^111^121) ^ . (2.4) 

V ^12 S12 - Sj2 J 

Inserting the factorized expression for the phase space measure in (2.1), the subtrac- 
tion terms can be integrated over the antenna phase space. In the case of initial-initial 
subtraction terms, the antenna phase space is trivial: the two remaining Dirac delta func- 
tions can be combined with the one particle phase space, such that there are no integrals 
left. We define the initial-initial integrated antenna functions as follows: 

Xik,j{xi,X2) = YTTT / [dkj]xiX2d{xi -xi)6{x2 -X2)Xikj, (2.5) 

where we introduced C(e) = (47r)^/(87r^)e~'''-s^. 

Substituting the one-particle phase space, and carrying out the integrations over the Dirac 

delta functions, we have. 



,2,-e -^^" 



Xikj{xi,X2) = {Q ) ^ . _ . J{XI,X2)Q XikJ , 

with Q^ = q^ = {pi + P2 — kj)"^. The Jacobian factor, JT{xi,X2) is given by 



(2.6) 



^f N XlX2(l+XlX2) ^^ f {l + Xi){l+X2) \ ' 

and the two-particle invariants are given by: 

_ xi (1 - xl) _ X2 (1 - xl) 

Sij — —Si2 ; , Sj2 — —S12 ■ . (,^.OJ 

Xl + X2 Xl + X2 

The integrated subtraction term is then, 

^,s,n, = ^^_^ /^^C(e)^..,(x„X2) 

^, ^-^ Jm+l J Xl X2 
m+l J 

Xd^rn{ki, ..., kj^i, kj+i, ..., km+i;XiPi,X2P2) 



xjj^\ki, . . . , kj^i, kj+i, . . . , km+i) , (2.9) 

where we have relabeled all ki —^ ki. The final step is to convolute this subtraction term 
with the parton distribution functions of the initial state particles. The integrated version 
of the subtraction pieces is then combined with the virtual and mass factorization terms 
to yield a finite contribution when e ^- 0. Recasting the convolutions appropriately, the 
integrated subtraction term is 

^^^S^+lJ ii J i2 7gi Xi 7^2 X2 \X2j \X2 

xC{e)X,k,j{xi,X2)<laB{iiH^,^2H2). (2.10) 

This convolution already has the appropriate structure and combination with the virtual 
corrections and mass factorization can be carried out explicitly leaving a finite contribution. 
The remaining phase space integration, implicit in the Born cross section, do"^, and the 
convolutions can be safely evaluated numerically. 

At NNLO, two types of contributions to m-jet observables require subtraction: the 
tree-level {ra + 2)-parton matrix elements (where one or two partons can become unre- 
solved), and the one-loop (m -|- l)-parton matrix elements (where one parton can become 
unresolved). The corresponding subtraction terms are denoted by da^^j^^ and da^y^^^Q. 
Antenna subtraction terms for the final-final [27] and initial-final [46] cases have been de- 
rived previously. In the initial-initial case, da^^j^Q was derived in [45,47]. It contains 
subtraction terms for single unresolved limits (each containing a single three-particle an- 
tenna function X^) and for double unresolved limits (containing four particle antenna 
functions X^, products of three-particle antenna functions {X^ ■ X^) and soft large-angle 
correction terms {S ■ -^^3)). 

The integrand in the (m -|- l)-parton channel consists, besides the one-loop (m + 1)- 
parton matrix elements, of several contributions (independent of whether the radiators are 
in the initial or final state): 

(a) The integrated one-particle unresolved subtraction terms from the [m + 2)-parton 
channel, which cancel the explicit infrared poles of the virtual one-loop (m + 1)- 
parton matrix element. 

(b) The virtual-unresolved subtraction term da^^^Q which subtracts all single unre- 
solved limits from the virtual one-loop [m + l)-parton matrix element. 

(c) Terms common to both above contributions, which are over subtracted. Each of these 
terms is formed by a product of an integrated and an unintegrated three-parton tree- 
level antenna function [X^ '^z)- These terms contain the full set of singly integrated 
(X3 •X3)-terms from da^jy^^Q, plus additional terms which must be further integrated 
down to the m-parton channel. 



(d) The integrated soft large-angle correction terms {S ■ X^). 



(e) Terms arising from the mass factorization of the parton distribution functions at 
NLO. 

Unintegrated subtraction terms newly introduced in the {m + l)-parton channel have to 
be compensated by their integrated forms in the m-parton channel. The integration of 
contributions of type (b) in initial-initial kinematics is the main topic of this paper, they 
are already known for final-final [27] and initial-final [46] kinematics. It should be noted 
that integration of terms of type (c) does not require any new integrals beyond the X^ 
already needed at NLO [44]. In particular, those terms obtained by integrating (Xg • X^) 
from dcj^jy^Q depend on the full set momenta of the (m -|- 1) partons in a non-factorizable 
way, but are not integrated any further. Any additional terms of type (c) are chosen such 
that the X^ depends only on m-parton momenta obtained from the phase space mapping, 
such that the integration of {X^ ■ X^) factorizes, involving only the known integral of X^. 
Terms of type (d) will be dealt with elsewhere [53]. 

With radiator partons i and k in the initial state, the contribution of type (b) reads: 
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(2.11) 



In here, X^^ ■ denotes a one-loop three-parton initial-initial antenna function, which is 
the only new ingredient. These antenna functions can be obtained by crossing from their 
final-final counterparts, listed in [27], and have to be integrated over the appropriate phase 
space according to (2.6). 

3. Integration of one-loop antenna functions 

The one-loop antenna functions are derived from one-loop squared matrix elements for 
all 2 ^- 2 processes [28] obtained from 7* — ?• qqg (quark-antiquark antenna functions), 
X ~^ 999 and x ~^ 9QQ (quark-gluon antenna functions) and H — > ggg, H — > gqq (gluon- 
gluon antenna functions) by crossing the off-shell current into the final state and two 
partons into the initial state. We denote a three-particle initial-initial antenna function 
with partons {i,j) in the initial state and parton k in the final state as X^- f^ at tree-level 
and as X^- ^ at one- loop. The tree- level and one- loop initial- initial antenna functions are 
summarized in Tables 1-3. The usual notation is used ,i.e. X^- ^ for the leading colour (iV) 
term, Xj- ^ for the subleading (1/A^) term and Xj- ^ for the Np part. 



Quark-antiquark initiated 


tree level 


one-loop 


quark-quark 
qq^ g 


40 
qq,g 


4I 4I 41 
qq,g^ qq,g' qq,g 



quark-gluon 

1Q ^ Q ^qq',q' ^ qq' ,q' ^ ^ qq' ,q' ^ ^ qq' ,q' 

1 ^ ^ 1 ^q'q',q ^q'q',q' ^q'q',q' ^q'q',q 

gluon-gluon 

yy y qq^g qq^g' qq^g^ qq^g 

Table 1: List of tree level and one loop three- parton antenna functions for the configurations with 

a quark-antiquark system in the initial state. 



Quark-gluon initiated 


tree level 


one-loop 


quark-quark 






qg^q 


40 


41 4I 41 
qg,q^ qg,q' qg,q 


quark-gluon 






Q9-^ g 


^qg,g 


qg,g^ qg,g 


gluon-gluon 






qg^q 


q9,q 


(^1 ril ril 

qg,q^ qg,q qg,q 



Table 2: List of tree level and one loop thrcc-parton antenna functions for the configurations with 
a quark-gluon system in the initial state. 



Gluon-gluon initiated 


tree level 


one-loop 


quark-gluon 
99 ^q 


^gg,q 


^gg,q' ^gg,q 


gluon-gluon 
99^ 9 


gg,g 


7?1 pi 
99.9' gg,g 



Table 3: List of tree level and one loop three-parton antenna functions for the configurations with 
a gluon-gluon system in the initial state. 



We start from the unrenormalized one-loop squared three-parton matrix elements (nor- 
malized to the corresponding two-parton matrix element and divided by a normalization 
factor C(e)) relevant to a particular antenna function, which we denote as X^-^^. The 
antenna function is obtained after renormalization and subtraction of the corresponding 



tree-level antenna function multiplied by the one-loop correction to the hard radiator pair. 
Renormalization of the one-loop antenna functions is always carried out in the MS-schenie 
at fixed renormalization scale /i'^ = Q^. It amounts to a renormalization of the strong cou- 
pling constant and (in the case of the quark-gluon and gluon-gluon antenna functions) to a 
renormalization of the effective operators used to couple an external current to the partonic 
radiators. The relation between renormalized and unrenormalized one- loop squared matrix 
elements is as follows: 

^i,R _ ^i,u bo Vo yo .on 

^ij,k - ^ij,k ~ ~ ^ij,k - — ^ij,k ) 1-3- -Lj 

4'l = 4'S ' (3.2) 

y1,r _ y1,u bo,F yo 'nojn yO ,„ „x 

where 

^0 = ^ , bo,F = -l (3.4) 

are the colour-ordered coefficients of the one-loop QCD /3-function: 

/3o = boN + bo^pNF . (3.5) 

The renormalization constants for the effective operators are 

110 = 0, ?7o,F = ior X = A, 
3 
■no = bo + -, rio,F = &o,F ior X = D,E , 

r?o = 2 6o , ilo,F = 2 6o,F for X = F,G. 

The one-loop antenna functions are obtained from the renormalized one-loop squared 
matrix elements by subtracting from them the product of the tree-level antenna function 
with the virtual one-loop hard radiator vertex correction [6,27]: 

^ii,k = ^ii.k ~ '^2 ^ij± ' (3-6) 



-^ij,k — ■^ij,k ~ ^^2 ^ij,k 1 (3-7) 

^ii,k = ^ij.k ~ '^2 ^i,jk ■ (3-8) 



The one- loop corrections to the hard radiator vertex are listed in [27] for A}^, T>\, 7)^, J-'g 
and ^2- From these, the remaining functions follow: 

Al = Al, Al = Q, (3.9) 

P^ = 0, (3.10) 

8l = Vl 41 = 0, £l=Vl, (3.11) 

7-1 = 0, (3.12) 

^2^=7-2^ Ql = Q, Q\=fl. (3.13) 

The integrated forms of the single real radiation antenna functions are still differential 

in x\ and X2, such that no explicit integration has to be carried out. However, endpoint 



singularities can occur in either or both of these variables, which have to be regularized 
dimensionally. This regularization is obtained from the d-dimensional integrated antenna 
functions (2.6) by expanding the product of the Jacobian factor and the antenna function 
in distributions using 



(1 



-l-e 



.l,a 



z + 



n 



-VJl 



n\ 



(3.14) 



with 



P„(l - z) 



In" (1 - z) 
1-z 



The expansion is straightforward for the tree- level antenna functions [44] , which contain 
only rational factors in the invariants 



1 



Sl2 



X\X2 



'li 



X2(xi +X2) ' 



S2j 



x\ 



Xi(xi +X2) ' 



(3.15) 



where we denote the pair of initial state partons with indices (1,2) and we refer to the 
unresolved one with the index j. The one-loop antenna functions contain logarithms and 
polylogarithms, yielding branch-cuts at the kinematical endpoints, which forbid a direct 
expansion in distributions. Instead, we start from the unintegrated expression in terms of 
one-loop master integrals. Only two types of master integrals appear: the one-loop bubble 

"(47r)^r(l + e)r2(l-e) 



e(l-2e) 



Bub(n , — „ „, 

^^ ' 167r2 r(l - 2e) 

and the general one-loop box with one off-shell leg 

2(1 



-A 



A 



2,LO 



-P 



2\-^ 



(3.16) 



Box(sij,SjA: 



26), 1 



+ 






2F1 



2F1 



-e;l 



e; 



^ijk 



->%] 



Sik 



^ijk 



Sik Sijfi; ^ 

Sijk^ijSik 
\Sij Sijl^jySik Sijk) 



-e;l -e; 



2F1 



Sik 



Sijk Sik 



-e;l - e; 



SijkySijk 



Sik) 



(Sijk SijjySijjs: SUf^j 



(3.17) 



Both integrals appear in all kinematical crossings. The expansion of the terms involving the 
bubble integral is trivial, while the one-loop box Box(sjj, Sj^) contains the rational factor 
l/(sij Sik) a-nd appears in the unintegrated antenna functions with further rational prefac- 
tors. In terms of the expansion in distributions, one has to distinguish three prefactors: 1, 

Sij/Sjk and Sik/Sjk. 

In order to analyse the initial-initial kinematical configuration we limit ourselves to 
the case 



Sl2> q > 0, Sij < 0, S2j < 0, 



(3.18) 



and we study all the possible crossings of the master integrals. For each of the three 
crossings of the box integral, expansions have to be derived for each of the three prefactors 
mentioned above. 

These expansions proceed by analytic continuation of the hypergeometric functions in 
the box master integrals to the appropriate region of analyticity, and requiring that the 
limit Xi ^^ 1 does not result in an argument of the hypergeometric function equal to 1 or 
infinity (avoiding the branch-cut). In this situation, (3.14) can be applied safely to the 
coefficients of the hypergeometric function. 

Using the following notation 

Mbox (sii,Sjfc,-^ I =y^i-rJ{xi,X2)-^TZ(Box{sij,Sik)), (3.19) 

V Spq J Ly (^ej Spq 

Mbub(p') =-^J{xi,X2)n (Bub(p2)) , (3.20) 

where TZ selects the real part, we list the master integrals with the relevant prefactors in 
Appendix A. 

The resulting expressions for the integrated antenna functions X^- ^ are very lengthy, 
such that we only quote one of them, T)qg^g, in Appendix B as an example. Analytic 
expressions for all of them, as well as for the tree-level antenna functions Xj^- ^ expanded 
through to O(e^) are attached with the arXiv-submission of this article. 

4. Conclusions 

In this paper, we extended the antenna subtraction formalism to handle single unresolved 
radiation at one loop for processes with two partons in the initial state, as required for 
hadron collider cross sections at NNLO accuracy. The corresponding virtual unresolved 
subtraction terms consist of tree-level and one-loop antenna functions with both radiators 
in the initial state. These initial-initial antenna functions are required in unintegrated and 
integrated form. The unintegrated initial-initial antenna functions are obtained straight- 
forwardly from analytic continuation of the corresponding final-final antennae. The inte- 
gration of the one-loop antennae over the phase space relevant to the initial-initial configu- 
rations requires an expansion in distributions around the kinematical endpoints of the two 
initial-state momentum fractions, which we performed for all relevant master integrals. 

Using the results of this paper in combination with the one-loop antenna functions in 
the initial-final [46] and final-final [27] case, the NNLO subtraction terms for the one-loop 
{n + l)-parton contribution to n-jet observables at hadron colliders can be constructed 
and implemented. To accomplish a full NNLO description of n-jet observables at hadron 
colliders, subtraction terms for the double real radiation (n + 2)-parton contribution are 
equally needed. The construction of these subtraction terms in the case of hadronic colli- 
sions has been described in [45]. A large fraction of the integrated antenna functions have 
been derived already [47]. Once the full set of integrated double real radiation antenna 
functions is completed for the initial-initial case, the antenna subtraction method can be 
applied to the computation of NNLO corrections to jet production at hadron colliders. 
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A. Master integrals in initital-initial configuration 

The initial-initial one-loop antenna functions can be expressed by the bubble and box mas- 
ter integrals defined in (3.19-3.20). They require analytic continuation to the appropriate 
kinematical region (3.18). 

The analytic continuation of the general bubble integral is straightforward 

Mbub(p') = ^^J{xuX2)A2,Lon{-p^ -i5)-\ (A.l) 

for p^ > 0. The box integrals read: 



57r2 ^,, _,, .1 



(g2)2+^Mbox(si„ S2„ 1) = -^6(1 - Xi)6il - X2) - -^<5(1 - X2){1 + Xi - Po(l - xi))- 



96e2"v- -"./"v^ ■".; - 2e2 



b{\ - xr)b(\ - xa) + -^ (2 + 2x1 - 2?o(l - x^))Vq{\ - x^)- 



8{\ - X2) Uv^{\ - xi) - 2(1 + xi) log(l - xi) - ^^ log Qa;i(l + ^^i)') )) 
((1 + 2:2)^ + 2x1(1 + X2)2 + x\{\ + 2x2(1 + X2))) - 



2e2(l + a;i)(l + X2) 
1 



2^ 24(2(1 + xi)2?i(l-X2)+ Po(l- 2:2) ( -2Pi(l-xi)+ 

2(1 + xi) log(l - xi) + ^^ log Qxi(l + xi)2 j j j + ,5(1 - X2) (^ - 5^2(1 + ^^) + 

5^22?q(i _ x^) _ 24P2(1 - 2:1) + -^— ( - 4x2 iQg(4(^ _ ^^^^ ^^^^^ _ ^^^^ 

1 -xi V 

41og2(l - xi) + 4x? log(l - xi) log (^(1 + xi)2^ + X? (^log2 (^^^^-^^ - 2(^log2 Q(l + xi)' 

^ -((l+.X2)2+2xi(l+X2)2+x2(l + 2x2(1 +X2)))log(l-Xi) + 



e V(l + 2;i)(l + .^2) 

X\ ( A (^ , ^2^ 2.X?x|(l+XlX2) /xi(l+XlX2) 

— -log 4xil + xi -7-^ 2T71 2Tlog T 

(1-Xl)(l-X2) V / (1 -2;f)(l -X^) V 2:1 +X2 

(X1+X2) , , 4X1X2(1 +X1X2) /(2;i+X2)2 



log(4)-^F:^^^7T^^^log 



(1-Xi)(l-X2) ^' ' (l-x2)(l-x2) '^V (1+2:1) 

2xfx|(l + xiX2) , / (l+2;i)2 



1 



(l-x2)(l-x2) '' V (2^1+ 2;2)(1+ 2:1X2), 

ViTx^5(\ - xi)(5(l - X2) - 160 ( - 48 1 Pi(l - xi)2?i(l - X2)- 



3840 
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(1 + a;2)2?2(l - xi) -Viil- X2) [ 2(1 + x^) log(l - xi) + j^ log Qxi(l + xif ]]] + 
Vo{l-X2)[5n'Vo{l-Xi) + 2[ - 57r2(l + xi) + 24-°^ (1 - ^^i) , 



1 -X 



6[ -4V2{l-x^) + -^(log^(—^-r]+4\og{l-x,)\or^^^^'''^^ 



2(V('l(l + .i)^')+log(-i)log('^^,i^')V2Li2(l-:.i)mV 



4' '' J ^^ " ^ V4(l-a;i), 

(5(1 - X2) { 107r2X'i(l - xx) - 161)3(1 - xx) - IOtt^ log(l - a;i)(l + a;i)+ 

1 



161og3(l - .Ti)(l + xx) + —^ - 7r2 log(a;i) + 241og2(l - xi) log 

1 - .Tl V 

n /(l+.Tl)2\ ^ X 9 / 4 \ o 

_■: 1 I *- ^^ 1 I 101 /I ^ \1„„^ I 1 I 01 i 



I-.T1V ^- ' ^ ' ^V4(l + 2;i) 



^2 



5.^ log ^^-^ j + 12 log(l - .,) log^ ^ (^^^ j + 2 log ^ ^^ ^ ^^^ 

241og(l - xi) log2 (J(l + xif\ + 41og3 Q(l + a;i)2 ) + 481og2(l - x^) log ( (1 + x^f ) + 



28(1 + a;i - Po(l - xi))C3 - 777^ ^, 7 log 



(l+a;i)(l+a;2) 



2(l-a;i)(l-a;2) " \{l + xi)\ 
((1 + X2f + 2x1(1 + X2Y + xl{\ + 2x2(1 + X2))) log2(l - Xl)+ 



l0g(xi)l0g ( ^^jP^^ +2 log f ,, ,^^ ,, ) (-log(l - Xl) - log(l - X2)) ) + 



(l-xi)(l-X2) V ^V4(l-a;i); ^ \{l + x^) 

8x^x|(l + xiX2)log(l - xi) / (X1+X2)'' \ x^X2(l + xiX2) 2/ (a;i+X2) 



log TT^-^TTT^^^-T^ - 7T^^l2T7r^-^log 



(l-x2)(l-x2) ^V(l+ 2^1)2(1 +X2)V (l--^?)(l-2;i) "" V(l+2;i)'(l+^2)2 

((1 + X2f + 2xi(l + X2)2 + x2(l + 2X2(1 + X2))) log(l - X2) l0g(l - Xi) + 



(1+Xi)(l+X2) 

- 2x^ log(x2) + 4((x? + xl) log(4) - 2x1 log(l + a^i) - 2x^ log(l + X2))+ 



log(l-Xi) / ^ 2,„ . . , ,,. 2 , 2x,„„,,x o^2,„„,i,^x o^2- 



(1-Xi)(l-X2) 

, x?xi(l + X1X2) ^^g / a±xini±^\ \ ^ /4xfxi(l + xiX2) j^g^^ _ ^^^ 



(1+Xl)(l+X2) \{xi+X2Y{l+XlX2f J J \{l ~ X{){1 - Xi) 

2x^X2(1 +X1X2) / X2(xi+X2)^ \\ /xi(l+XiX2) 



(1 - x2)(l - xl) ^°^ \{l- xl){l + X2Y JJ^^'^y Xy+X2 
2A(I^X^)(\^X') K(^(^ " ^2)(1 + 2;i)(l + X2)2 + X?(-5 - 5X2 + 6x3) + 
x\{-h + X2(-5 + 6X2))) + 12(1 + Xi)(l + X2)(x2 + xl) log'(4) + 



„2,_/ l+^Ml„„.i , . N , .2, (^^X2 



48(1 + Xi)(l + X2) Xf log —-^ log(l + Xl) + X^ log —-^ l0g(l + X2) + 



„2 2,, , ^ ^ x,„„2^ (1+Xi)2(l+X2)^ 



12xfx^(l + X1X2) log" -^ ^^-^ ^ + 12xf (1 + Xl)(l + X2)Li2(l - Xl) + 

V(xi +X2)'^(1 +XiX2)'^/ 

12x2 ((1 _^ ,^^^(^ _^ X2)Li2(l - X2) + 2x2(1 + X1X2) ( Li2(-^'^^ ^ '^'^' ^ 



■(1 +XiX2)2 
^. X^_X|X2 . X^-X1X|\\\ 

Xl + X2 a;i + •^2 / / / 
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(g')'+^Mbox(si2, si„ 1) = - 2^a:2<5(l - X2) + ^xj ('22?o(l - 2^2) + S{1 - x^) log f^^—^\ 

- xl{x2 + 2x1 + 2:1 1 + 2X2 1 + X2+ xl))) ^\ i^ - ^-^1^ + 

eV (1 + a;2)(a;i +a;2) 1 - X2 

2xl{x2 + 2x1 + x,{l + 2X2(1 + X2+xl)))-^^^\^^^^--^( log ^^^±^' 

(1 + a:;2)(a;i +a;2) 1 - a;2 V 

xi(l+XiX2) / a:ia;|(xi+a:2)^ \\\ 

(xi + X2)(i + X2) °^ ui + ^i)(i + ^2)vy y ^ 

l^xl(^A8(v2il - X2) + Pi(l - X2)log (^^il-^^^ - 22?o(l - X2)(^^2 + 121og(xi)log (^^ 
6 lo~ ' 



3 log(x 
log 



(48(^2(1 -X2) + Pi(l-X2)log(^^|l^-^))-2Po(l-X2)(.^ + 121og(x01og(^}^^ 
-^^V-7r2 + 21og2(4) + 21og('l(l + xif')log('(l + xi)^)) + 



(1+Xi)2 

121og (^ll±^^Li2(l -xi) + 12Li3(l -xi) +28C3)) - ^^log(xi)log (^^)- 

X2(X2 + 2x2 + ^^(1 + 2X2(1 + X2 + xl)))—^^^^SL_IlL_ 

' 2(1+X2)(xi +X2) 

2X?(X2 + 2x2 + Xi(l + 2X2(1 + X2 + x2)))-i^lii^^H) + log(l - Xi) f 4 

' (1 +X2)(X1 +X2) V 1 



log(xi) 



(X2 + 2x^ + 



2-2 / 

7 ■ ^7^ 2T (1 + ''2){xi + X2) log(4) - 2(1 + X2)(xi + X2) log(l + Xi)- 

(Xl +X2)(1-X^) V 

^ S/-, \ , /^ XiX2(xi+X2)^ \\\ , ,_, ,/2x?l0g(xi) ^ 2/ 

2x3(1 + X,X2) log ( (TTibaT^J J J + ^'^^^^ " "^^ l^^ - '"^^^ - 

Xi(l + 2X2(1 + X2 + X^))) g\ ]> - ' (1 + X2)(X1 + X2) log(4)- 

' (1+X2)(xi +X2) (Xi +X2)(1 -X^j V 

2(1 + X2)(X1 +X2)l0g(l+Xi) + 2x^(1 +xiX2)log(^^^l^|^i±^)))- 
12(x,+X2)(l-xl) "K-"'(^-"^H^2 + 2x2+x,(l + 2x2(l + X2+x2)))- 
6(1 + X2)(xi + X2) log2(4) + 24(1 + X2)(xi + X2) log (^^^ ) log(l + xi)- 

TO 3n , \f^ 2/ X2(xi+X2)3 \ 2 /2;2(X1 +X2)3(1 -X1X2)' 

12X9(1 +X1X2) log ^ — —-rr + log — 

'^ \ ^ \{1+Xi){l+X2)y ^ \ (1+Xi)(l+X2)2 ^ 

/ X9('xi + X-^)^ \ /xiH +XiX9'l\ r, / (xi + X9^2f1 — x?X^ 



^]nrf ^2(a;i+X2)3 \ / xi(l + XiX2) \ 9 f {xi + X2)Hl ~ xlx?.) \\ 

^{{i+x,){i+x2r) H 2;i+x2 ) ^ \ (i+xi)(i+x2)2 ;; 

12(1 + X2)(X1 + X2)Li2(l - Xi) + 24x^(1 + X1X2) ( - Li2 (-, ,^^~,^'^'^' r ) + 

V V (2^1 +2;2)(i -2:1x2)/ 



13 



Lis 



X2 — x1x2^ 


)+Li. 


^(1- 


-xf)xi 


a;i + X2 y 


/ 


V i- 


~ 3^1 Xq 



Oie) 



(A.3) 



Sl2 



(Q^)^+^Mbox(si2, Si,, — ) = -TTlH^ - ^i)'^(l - ^2) - TTl^il - 2;2)(1 + xi- 2?o(l - a^i)) 



S2j 



2e 



2e3 



^<5(1 - a;i)(-l - X2 + Vail - X2)) + ^2?o(l - a;i)(l + X2 - Po(l - 2:2)) + ^(1 + xi)Po(l - X2) 

gO g^ gZ 

^(5(1 - X2) ( - X>i(l - xi) + —^ + -; — i— log 

1 



1 — Xi 



1 — a;! 



4(1 -xi) 



24e2 



(5(1 -xi) - 7r^(5(l - X2) - 24 -22?i(l-a;2) + 2(l + .T2)log(l-a;2) + 



^2 i„„ /' 1 2^ 



log -X^(l+X2) 

1 - a;2 V 2 



24e 



1 (1 + X2)^ + 2X1(1 + X2)^ +xl{l + 2X2(1 + X2)) 



(1+Xi)(l+X2) 



- 24 - (l + X2-Po(l-a;2))I?i(l-xi 



2(1 + xi - Po(l - a:i))Pi(l - X2) + Po(l - 2^2) ( (1 + ^i) log(l - xi)+ 

^^ log Qxi(l + xif))) - 2(1 + .T2)2?o(l - xi) log(l - X2)- 
■^2 Ti_^i „M^„ / ^„2n I „ ^ M xn „ ^/ „2/ 



^_^^-2?o(l-a;i)log(^-x^(l + X2)) ) -6il-X2)[ -n'{l + x,)+ 



TT^Voil - xi) + 6 



1 — Xl 



+ 6 -P2(l-.xi) 



1 — Xl 



log' (4)+ 



21og(4)log 



(^^) +^°^'^^ " ""'^ " ^^°^^^ " ^i)log(a;i) - 2 log (^4(^1-^^^ log (^ 



(l + a;i)M + 



log2 ( (1 + a^i)M - 2Li2(l - Xl) j j j j + 2<5(1 - xi) ( - 7r2(l + X2) + 7t^Vo{1 - X2) 
24P2(1 - X2) + 24i^^4ii-^ - 6^^^ ('x2 ('41og(l - X2) log f ^i^l^^ 



1 - X2 



1 - X2 



log' ( ^2^2(1 + 2^2)] + l0g(x2) log (^ix2(l - xlf 



1 



((1 + X2f + 2xi(l + X2)' + x2(l + 2x2(1 + X2)))- 



X2(l + X2) 
US{l-X2)C3]) + 
1 



■log(l-xi)+ 



(1+Xi)(l+X2) 
2((1 + X2)' + 2xi(l + X2)' + x2(l + 2X2(1 + 2:2)))-— -7— r log(l - X2) + 



(1+Xi)(l+X2) 
Qx2(l+X; 
X2(xi+X2)'^ \ 2x^X2(1 +X1X2) /xi(l+XiX2) 



(1-X1K1-X2) (^^ ^°s [h^' + ^^^V + ^' ^°s (^'^^ + ^') 

X^X2(l +X1X2) 

2— ^tt:^ ^loa 



(l-x2)(l-x2) ^V(l+2;i)(l+ 2:2)2(1 +X1X2); (l-x2)(l-x2) 



■log 



Xl + X2 



— (12(1 + X2)2?2(l - 2:1) - 12(4Pi(l - xi)2?i(l - X2) + Vail - X2)p2(l - 2:1))+ 

48(1 + xi)2?2(l - 2:2) + Sil - X2) (iT^Viil - Xl) - 22?3(1 - x^)) + 48I?i(l - X2) log(l - xi)+ 

48x2 /I \ 

48xi2?i(l - X2) log(l - Xl) + ^I?i(l - X2) log -xi(l + Xl)' + 48Pi(l - xi) log(l - X2) 

1-xi V4 J 
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2Ax^ /I \ 2 / 

48a;2Pi(l - xi) log(l - X2) + 2_p^(i _ xi) log -a;^(l + X2) ] - Po(l - xi) 7r2(l - xj)- 

l-a;2 V2 /I- 2:2 V 

24(1 - xj) log2(l - X2) - 24x2 i^g(^ _ ^^) i^g / 1^2(^ ^ ^ ^ g^2 / i^g2 / 1 ^^(^ ^ ^^^ j ^ 
log(x2) log Qx2(l - ^2)'))) + Y^^o(l - '^2) (" - 7r2(l - x?) + 6(1 - a;2) log2(l - xi)+ 
12x1 log(l - a;i) log(a;i) + &xl log ' ^^^ " ^'^ 

-^25(1 - xi) f - 16(1 - xl) V(l - X2) - 24x2 iQg2^^ _ ^^^ j^g j" 1 2(^ ^ ^^) 
I-X2 V V2 

2 log(l - X2) L^{1 -xl)+ 6x2 l" i^g2 Q^^(^ ^ ^^)A ^ i^g(^^) i^g Q^^(^ _ ^2)2\ \ V 

x2 (^^2 log Q^^(i + ^^)^ _ 2 l0g3 Qx2(l + X2)^ - 77r2 log ^(1 - X^)^ + 2 log3 Q(l - x2)^ + 

7^2 log Qa;2(l _ xl)\ - 2 log3 Qx2(l - xl)\ \ - 28(1 - xDC^ + 22?o(l - a^i) [j^^Va{l - xa)- 



.?(! + . 0^/°Hi^' + ^^)')+^'^?^^^(^-^^: 



24I?2(1-X2)- 14(5(1 -X2)C: 



1 



5(1 - X2) 2x2 iog^(4) ^ iQg(^ - Xl) U 



1 — Xl 

2 log(l - xi)(x2 log(64) - (1 - x\) log(l - Xl)) - 6x2 iog(i6(l - xi)) log(xi) ) + 



6xi log(4) log 



1 — Xl 
Xl 



log( 



— --4 ) + Xf TT^ - 6l0g^(4) - 6 log(l - Xl) log , 2 

yi + Xl j \ \ Xt 



1 — Xl 



X log ( (1 + xi)2 ) - 2x? log3 ((1 + x^fj - A log (4^^) (7r2 - 61og2((l + xi)2))- 
28C3 + 4x2 I 3 log I il±^V i2(i _ xx) + 3Li3(l - Xl) + 7(3") \ + y^<5(1 - Xi)(477r45(l - X2)+ 



480(^2p^(l _ ^2) _ 8Pg(l _ ^2) _ 14P^(1 _ ^2)^_^)) ^ 



2x2 



(1-Xi)(l-X2) 



log (xi) log 



1+Xi 



((1 + X^f + 2xi(l + X2)2 + X?(l + 2X2(1 + X2))) 



log2(l -Xl) 
2(1+Xi)(l+X2)" 



2((1 + X2)'' + 2xi(l + X2)^ + x2(l + 2x2(1 + X2))) 



\0i{l-X2) 
(1+Xi)(l+X2) 



^^^^^y^^^i^ (^(1 + Xi)(l + X2) (^X? log Qxi(l + Xi)2^ + x2 log Qx2(l + X2) 



2X1X2(1 +XiX2)log 



x|(xi + X2)'' 



(l+Xi)(l+X2)2(l+XiX2) 

2 I O/n (a I ri^ \2 I ^2 



l0g(l-X2)( -2((l+X2)2+2xi(l+X2)2+x2(l+ 2X2(1 +X2)) V/°g^\.^';^^ . - 

V V y (i + xi)(i + x2) 

,-^ _ 2,-^ _ ^2) ((1 + 2;i)(l + X2) [x\ log (^ixi(l + Xi)2 j + x2 log Qx2(l + X2) 



2X1X2(1 + X1X2) log 



X^Xl +X2)'' 



/ 2X1X2(1 +XlX2)log(l - Xl) 



Xl+a;i)(l+X2)2(l + xiX2);;; ' V {"^-x\){\-xl) 

4X^X^(1 + X1X2) log(l - X2) 2x\x\{\ + X1X2) / X2(X1+X2)'^ \\ /^Xl(l+XlX2) 



(l-x2)(l-x2) 



(l-x2)(l-x2) ^\{\^X^){\+X2) 



Xl + X2 
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1 



TT^il - Xl)(l - X2) (1 + X2)'^ + 2xiil + X2f +xl{l + 2X2(1 + X2)) + 



I2{l-xl){l-xl) 

6(1 + a;i)(l + X2) i^xl + xl ] \og^{2) - 48x^(1 + xi)(l + X2) log(2) log(l + a;i) + 



24xf (1 + a;i)(l + X2) log"(l + xi) + 6x^ - 2(1 + xi)(l + X2) log(2) log(x2(l + X2)) + 



(1 + Xi)(l + X2) l0g"(x2(l + X2)) - 2xf (1 + X1X2) log 



2 ( X2(xi+X2)^ 
(l+.Tl)(l+.T2)2 



+ 



log' 



X2(xi +a;2)^(l - X1X2) 



(l + x,)(l + x.)^ ^^-(l + x,)(l + ...)log^^2(^-x^)) + 



1 



2xi(l + xia;2)log 



(1 -XI)X2 



I2xi\ - (l+Xi)(l+X2)Li2(l -Xi)+ 2x^(1+ X1X2) -Li2 



(Xi +X2)(1 -X1X2) 



^. X2-XIX2X ^.J-(1-X^)X2 



Xl + X2 



-L Uj 1 Xo 



0{e), 



(A.4) 



2x^X2 \ 1 



(Q2)2+-Mbox(si2, Slj, ^) = -^(5(1 - Xl) + :^ ( 2?o(l - Xl) - (5(1 - Xl) log 



^ -2 



1 + Xl Xl + X2/ 12e 
xn \l 2 R, 2/^2^2(1 + 2:2) 

0(1 — Xl) TT — 6 log ' 



-— x^ -12 Pi(l-xi)+Po(l-2;i)log 



f xli^+X2) 
l2(l-X2)2 

2(1 -X2)^ 
X^l +X2) 



Rl /Ml 1^2:2(1 +X2)2 

6 log ( — 1 log ' 



2(1-X2)V \X2j ^V 4(1 -2:2)2 

-(x^(:^2+x.(l + 2xi + 2(l + xi+x^)x2)) ,^|°g^;.-f ^ -r^(^°g(2rTT4) 

ey ' (1 + xi)(xi +X2) 1-xiV v2(-l + 2;2)'^/ 



xf(l+XiX2) /^ XiX^(xi+X2)^ 



(l + xi)(xi +X2) V(l + 2:i)(l-2:i) 



2\2 



^x'i ( 6232(1 - Xl) - 1223i(l - xi) log 



x^(l + X2) 
2(1-X2)2 



Po(l-2;i) hr'-61og 



3:2(1 + 2:2) , , „, , - 1 , 

^^^-— ^)+61og(-)log^ ^^^_^^, 



1^^^JX2(1+X2)^ 

X2, 



<5(l-xi)( -TTT^logl Jj±^)+21og' 



4(1-X2)2 
Sf 1+2^2 \ , ^_2, „/2;2(l+X2) 



2(1 -X2) 



7tt^ log 



2(1 -X2) 



21og^ i'-^^T^) +-^log('|#^V21og^ f^^^ ) +28C3 



2(1~X2)V V 2(1 -2:2)2, 

2, ,N^ log2(l-Xi) 2x^X2(1 +X1X2) 



2(1 -X2) 

2(x2 + xi(i + 2x1 + 2(1 + Xl + x2)x2)) ^./"^ ';, 7 . + ,7^-^^2\7 T'\ 

^ ^ 2(1 + xi)(xi +X2) (1 - xj)(xi +X2) 

/x^(l + X2) 



^ /xi(l +XiX2)\ , / X2(xi+X2)^ 

X log ^ ■ log 



Xl + X2 



(l+Xi)(l-x2)2y l-xi 



■log(l -2:1) log 



V2(l-X2)2 



xf(l+XiX2) ^f XiX^(xi+X2)^ 



1 



(1+Xi)(xi+X2)""'' V(l+2;i)(l-2;2)2;; " 12(l-x2)(xi+X2) ' 
Xi(l + 2xi+2(l+Xi+x2)x2)) -6(1+Xi)(xi+X2)l0g2 (^i-t^V 



(^-7r2(l-Xi)(x2 + 



(1 + X1)(X1 + X2) l0g2 {^^^^±pl ) + 6( (1 + X1)(X1 + X2) lo; 



2-2x2 

2 / 2:2(1 +X2) 



2(-l + X2)^ 
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2x5*(l+a;ia;2) log -— — 5-- + log 



log^ f ^^^r^^l^'l? ) ) ) + 24.3(1 + .,.2) f - Li:- ^ (^ - ^^^^^ 



(i + .Ti)(i-.Ti)2 ;;; ^^ ^^v V(^i + 2;2)(i- 2:1x2) 

2:2 -a;?a;2\ , ^. /(I - xl)x^ 

', + i^l2 — 2~' 

Xi + X2 J \ I - X^X'2 



U4^^-^ ) +Li2( ^^^^^ ) ) ]+Oiel (A.5) 



(g2)2+^Mbox(si„ S2„ ^) = -:^.t2,5(1 - xi) + ^xl f 22?o(l - xi) + ,5(1 - xi) log ^ ^^^ ^'^ 



si2^" 263-2^'^- ..,, 2^2-2 V' -.; . ^v^ ■^^^^"n-^2(l + .T2)2 

1 X2 



e^ (1 + xi)(xi +X2) 



(X2 + Xl(l + 2X1 + 2(1 + Xl + x2)x2)) 



4(1-X2)2\\ , 1 2^^^ ^ ^, 2 



22?i(l-xi) + Po(l-xi)log \^ , "\, +--x^5(l-xi) 57r^ + 61og' 



2:2(1 + 2;2)V/ 24e 2 ^ ^^V V(l + 2^2)2 



.2 



1.1 , ^ ^2\ , 1 o i„„Y (1 - ^2)2 \,_/ (1 + 2:2)2 



12 log^ (^-(1 + X2)^ j + 12l0g (^^— ^ j log (^^^^-^ j - 12Li2(l - X2) j + 

1/9/ / . 9n n\ log((l — Xi)(l — X2)) X^ , / 4X2 

-['4{^2 + Ml + 2x,+2il + x,+xl)x,)) gVxi)(x;+X2) -T^^^HoriF 

4(1+Xi)x2 /l+X2\ 2X1X2(1 + X1X2) / XiX2(xi +X2)'' 

I2 "^S 1 o ) + 7i I2T7Z I „ \ ^^S ' 



\-x\ ^\ 2 ) (l-x2)(xi+X2) '"V(l + 2:i)'(l+a;2) 

^x^ (^48P2(1 - Xl) + 48Pi(l - Xl) log (^^1-^^ + 2Po(l - 2:1) (^ - 5^^ - 6 \oi ^ ^^ ^ ^^^, 

4(l-2)^,.._2A. ,..2^,..l,g(,^)l,g^Il±M^ 

V4(l-X2 

4 \ „, ./ 4 



241og(l-X2)log( ^l^^j +121og^ (^-(1 + X2)2j +121og(x2)lo«^^^^_^^^ 



12Li2(l-X2))+,5(l-Xi)(l61og3(l-x2) + 57r2log(^^^--^j-21og3^^^^^^^^ 
10^2 log (1(1 + ^2)2^ _ 241og2(l - X2) log Q(l + X2)2') - 41og3 Q(l + X2)2 
21og(l - X2) ( Stt^ + 61og2 ( 7^^— 2) - 121og2 [-{\ + xa)^)) + log(x2) f^^- 



_2 



(1+X2)^ 
/(1+X2)2\\ ^^, /(l+2:2)2 



61og (77^ rr) ) -121og(^y^ ^)Li2(l-X2) + 12Li3(l-X2) + 28C3 

-2 ( 2(x3 + x3x2 + (-1 + xi)x3) log^(2)+ 



V4(i-X2)yy ='V4(i-x2) 

(x\ + xfx2 - (1 - xi)x3) log2(4) 

(1 -Xi)(xi +X2) ' (1 -Xi)(xi +X2)' 

xixi(l + X2) log ( i(l + xi)2 ) log((l - Xl)(l - X2)) + log(l - X2) ( - Xl(l - X2)X^ log(4)- 



4x3(1 + X2) l0g(l + Xl) + x3(l + X2) log U(l + Xi)2 j + x\ log Q(l + Xi)2 ) log((l - Xi)(l - X2)) + 

2:?X2 log \U\ + xi)2 j log((l - Xi)(l - X2)) - xix2 log { — {\ + Xa)" j - 

XixilogUl +X2)2 jj +xilog(l-xi)(^ -x2(l+Xi)(l+X2)l0gQ(l+Xi)2 ) + 

xl ("xi log(16) - 4(1 + xi)(l + X2) log(l + X2) + log U{\ + X2f\ + X2 log f 64(1 + X2)M - 
n \ y (a (^^ \2\\\\ 2x|x|(l+xiX2) /, ( 1 ^, ^ xi(l + xiX2) 

(l-Xl)x2l0g(^4x2(l+X2)jjjj- ^^_^2)(,^+,^) (^log(^(Y^^jlog(^ Xl + X2 



17 



\(l-Xiyj \ Xx+X-2. ) \ \Xx+x{X2) 

J_ 2/ 4 \ a:ilog(a:2) / (l + a;2)^ \ 2x3x1(1 + X1X2) 2/ (3^1+2:2)^ 

2-2x/°S V (1+2:2)2; 1-Xi '°^U(l-X2)y'^(l-x2)(xi+X2) ^^ ^ (1 + Xi)2(l + X2)2 

8x?x^(l + X1X2) log((l - xi)(l - ^2)) (i_,2)[,^^,^) log ( (TT^^)^(rf^) + 

x/ . 9/ /, ^ „/ 9n nn log(l — xi) 2x?(l + X2) log(xi) 

log 1 - X2 - 4x2 X2 + ^^ 1 + 23^1 + 2 1 + XI + X X2 ,, J'^ ., ]! , + -jf^ w^ Z V " 

V ' (l + xi)(xi +X2) (1 -xi)(xi +X2) 

■ ( - 4x3(1 + X2) log(l + xi) + x3(l + X2) log (4x1(1 + xi)2 ) + 

(1 -xi)(xi +X2) V V / 

xl [xi log(16) - 4(xi + X2) log(l + X2) + X2(xi log L^^^ j + log U(l + X2)2 j) n - 

X3x2 ^ (^+-^-2) / (1+X02(1+X2)2 \\ pXixi(l+X2)l0g(x2) _ 

^ ^{\-X^^{xr+X2) ^ V(a;i+ 2:2)2 (l+XiX2)2yy ^^ '\ (1-Xi)(xi+X2) 

(l-XiKx7lx2) (^?(^ + ^^) ^°^ ( (T^) + ^2 (-1 log(16) - 4X1(1 + X2) log(l + X2)- 
X2l0gf-X2(l+X2)2j +XiX2l0g f 4X2(1+ X2)2jjj - 4x^x2 —-L±g^^L_^ 
/ (1+Xi)2(l+X2)2 \\ 2x^x2(1 +X1X2) / X2(xi+X2)3 \ 

"" '°H(^i +^2)2(1 + xiX2)2y'y' + (1 - x2)(xi +X2) '°^ Ul - 2:?)(1 - AY) 

2xjx|(l +X1X2) / X2(l +XiX2) \ / Xi(xi+X2)3 \ 1 

(l-x2)(xi+X2) °^ ^ Xl + X2 j °^ lv(l-^?)'(l-a:i)j " 12(l-x2)(xi-, .,, 

X AW ( - 5x2 + Xl ( - 5(1 + X2) + Xi(-5 + 6x1(1 + X1X2))) j + 12(1 + xi)(xi + X2) log2(4) + 

48(1 + xi)(xi + X2) log (i±^) log(l + X2) + 12x3(1 + X1X2) log2 {,^^^^^^:) + 

12(1 + xi)(xi + X2)Li2(l - X2) + 24x3(1 + X1X2) f Li2 f i^l±^') - Li2 ("^l^^) _ 

-<^^))) --(')■ <") 

The remaining box integrals can be obtained from the previous ones by exploiting the symmetry with 
respect to the transformation xi o X2. More precisely, Mbox(si2, S2j, 1), Mbox(sij, S2j, 52^/^12), 
Mbox(si2,S2j, si2/sy) and Mbox(si2, S2j, S2j/sij) can be derived from Mbox(si2,si-,, 1), 
Mbox(sij, S2i, sij7si2), Mbox(si2, sij, S12/S2:,) and Mbox(si2, sij, Sij7s2j) respectively by exchanging 

Xl <-!> X2. 



B. Example of an integrated antenna function: ^L„ antenna 

The integrated one-loop antenna functions Xi^^j result in lengthy expressions. As an example of 
these, the function 2?^gg(xi,X2) reads: 



4xiX 
2 
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-l2\2£ Al ' - 1 „.. ^ ^.. ^ 1 



{QyD'^gJxi,X2) = ^6{1 - xi)S{l - X2) + ^<5(1 - a;2)(l + xi - 22?o(l - xi))+ 

^(5(1 - a;i)(-l + 2:2(2 + (-1 + .T2)a;2) - .T2Po(1 - X2)) + 

0X2 c 

-^ ( 6 { I (2x2(2 + xi) + xi (4 + .Ti(2 + a;i)(l + xD) X2 + 

36e \ V(l + 2^1)2^2(1 + a;2)(a;i +a;2)'^ V 

f 2 + a;?(3 + xi (10 + a;i(7 + 2a;i)) ) ja;^ + xi (9 + a;i(2 + xi) (9 + 2x1 (3 + xi))) x^ + 
f 3 + xi(l2 + xi (19 + 2x1 (7 + xi (5 + xi))) ) jx^ + M + xi(l + xi (4 + xi(2 + xi)) ) 

X 2(l + xi)x^ +2xi(3 + 2xi(3 + xi(2 + xi)))x^ + 2(l + 2xi(l + xi+x?))xZ ) - (2( h 

(2 - (1 - X2)x2))Po(l - xi) + (1 + xi - 22?o(l - a:i))I?o(l - 2:2)) + 6{l - X2) ("(1 - xi) + 
22?i(l - xi) - (1 + xi) log(l - xi) + ^1±^ log ( ^ 



1 — Xl \ 1 + X 

<5(1 - ^1) ( - ( - ^'-^(l - '^2) - 3(x2 - 4Pi(l - X2)) - 12(-^ + 

(2 - (1 - X2)X2)) log(l - X2)) - 12 ^^'/^'"^^"^^' log {-^ 

J (l-a;2)a;2 V 1 + 2^2 

1 ^0,2 », . ^^ n . ^ , 24(-l + X2(2 - (1 - X2)x2))Pi(l - xi) 



^^ ^ 27r2^(l - X2)I?o(l - Xl) + ^ ^ ^ ■^'■'" '' ^ + 12(1 + xi)2?i(l - X2)- 

72 V X2 

24(2?o(l - 3:2)^1(1 - Xl) + 2?o(l - a;i)2?i(l - X2)) - 125(1 - X2)2?2(l - 3:1)+ 

2<5(1 - Xl) U22?o(l - X2) - 62?2(1 - X2)\ - Y^^o(l - X2) ({l - xif + xl log(2(l - xi))+ 

6(log(4) + xl log(4(l - Xl))) log(l - Xl) - 61og2(l - xi) + 12 log(2(l - xi)) ({l - xi)2 + 

(l+x?)log(Y^))+61og(Y^)(2(l-xi)2 + (l + x?)log(Y^ 

12 / 

— Po(l - a;i) xl + 2(-l + X2(2 - (1 - X2)x2)) log(l - xa) 

X2 V 

24Po(l-xi) _ , ^^g / 2 \ _ 1 2^(1 _^J_ 6(1 _ ^^)(_i + ^^(2_ 

(1-X2)x2 VI + 2^2/ (l-a;2)a;2 V 

(1 - a;2)a;2)) log2(l - X2) - (-1 + X2) (6 + 2x2(-3 + 7x2) + n^i-l + X2(2 - (1 - X2)x2))- 

3x2 i^g j^_^_^^ + 6(1 - (1 - X2)X2)2 l0g2 (y^) + 3 l0g(l -X2){^-{1- .T2)x2 + 

4(l-(l-X2)x2)2log(YJY^))) + ^^^^,^^^^^;^^^^3^^_^,^ 12(x2(-2x2 + x3+ 



X^+X'jfl+Xa -2X^1 -Xf(l-X2)x2(-1+X2(l+X2)(2 + 5X2)) -Xl(l-X2)x^(5 + 

X2(6 + X2(4 + X2))) + x^( - 1 + X2 - 6x^ + x^(6 + log(8))) - x?( - 1 + X2(3 + xa + xj 
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&xl + 2a;^ - X2 log(8))) + a;f ( - 1 + 0:2(2 - lOa;^ + 3:2(9 + log(8)))) j + Kx2(l + 2:2 - 2x5)+ 

x5x2(l + a;2(4 + 3a;2 - 4a;^ - 4x^)) - (1 - X2)x^(2 + a;^(3 + 2(x2 + x^)))- 

x\{l - x^f (-2 + a;2(-3 + a;2(4 + a;2(7 + 2^2 + ^xl))) j + x^ - 4 + .T2 (6 + a;2( - 5+ 

a;2(-6 + .T2(2 + .T2 + ^x\))) j j + X1X2 ( - 4 + X2(6 + X2(-ll + x2(7 + 2x2(-2 + X2(2 + X2))))) j + 

X^ ( 2 + X2 ( - 3 + X2(4 + X2(8 - X2(7 + 2X2(-1 + X2 + 2x3)))) j j \ iQg(^ _ ^^)_^ 

X?X2f X^+xt(l+X2)+X?(l+X2)(l+3x2)+XiX^(3 + 4x2) + 3x?f X2+X^ j j log {- j + 

X2 f x-^(l + xi) + x?(3 + 4xi)x2 + 3x1 (1 + 2xi) (l + x?) x^ + (1 + 2xi)2x^ + (l + xi + xf) x\ j 

X log ( — ) + x^x^ log ( — ) + x?X2 {x\{\ + xi) + x-J(4 + X1)X2 + x^) log(2(l - X2))- 

( - 2x^x^ + x'^x3(3 - 4(x2 + x%)) - (1 - X2)x2(2 + x\{Z + 2(x2 + x^)))- 

x\(\ - xif f - 2 + X2( - 3 + X2(4 + X2(7 + 2x2 + 4x2))) J + a;2 / _ 4 _^ ^^^g ^ ^^(_5 _^ 2x2(-3+ 

X2 + 3x^))) j + X1X2 ( - 4 + X2(6 + X2( - 11 + x\{l + 2x2 (-2 + 3^2(2 + 3^2))))) ) + 

X^ 2 + X2 ( - 3 + X2(4 + X2(8 - X2(7 + 2X2(-1 + X2 + 2x3)))) j \ \ iQg(^ _ ^^)_^ 

2^X^(1 + X2) (1 + (-1 + X2)X2)^ + x3(l + X2)(l + (-1 + X2)X2)2 + XiX^(l + X2)(3 + X2)(l- 

(1 - a;2)x2)^ + 3x^x2 (l + x^)^ + x^ (l + X2 (2 + X2 (-2 + 4x2 - 2x3 + 3x^))) j log {-^ — — j + 

4x?x^ log f j - 2(1 + X1X2) f x^ + xjx^ + 

xiX2(2 + x^) + x^3(3 _ ^2 ^ 2x4) + x?(2 - x2 + 3x|) + xt(x| + 2x6)) log (^iii-^il^^^ 

^ 2-24 [2x\x\ + x?x3(-3 + 4(x2 + x%)) + (1 - X2)x2 (2 + x2(3 + 2(x2 + x^))) 

(1 -Xf)X2(xi +X2)'^(1 -X^) V 

x3(l - xif ( - 2 + X2( - 3 + X2(4 + X2(7 + 2x2 + 4x2))) j + ^^ U + X2( - 6 + X2(5 - 2x2(-3+ 
X2 + 3x^))) j + xiX2(4 - X2(6 + X2(-ll + x2(7 + 2x2(-2 + X2(2 + X2)))))) + xt ( - 2 + X2 (3+ 
X2( - 4 + X2( - 8 + X2(7 + 2x2(-l + X2 + 2x\))))\S\ log(l - X2) - 8(5(1 - xi)r5(l - x^)Q^ + 0(e). 



(B.l) 
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